Freezing into Stripe States in Two-Dimensional Ferromagnets and Crossing 

Probabilities in Critical Percolation 
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When a two-dimensional Ising ferromagnet is quenched from above the critical temperature to 
zero temperature, the system eventually converges to either a ground state (all spins aligned) or 
an infinitely long-lived metastable stripe state. By applying results from percolation theory, we 
analytically determine the probability to reach the stripe state as a function of the aspect ratio and 
the form of the boundary conditions. These predictions agree with simulation results. Our approach 
generally applies to coarsening dynamics of non-conserved scalar fields in two dimensions. 
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We investigate the fate of a kinetic two-dimensional 
Ising ferromagnet following a quench from above the crit- 
ical temperature to zero temperature. Intriguingly, this 
evolving system can sometimes get trapped in an in- 
finitely long-lived metastable stripe state that spans the 
system vertically or horizontally [H, [3, [^. The case of 
two dimensions is special; in one dimension the system 
always reaches the ground state, while in three or more 
dimensions the ground state is never reached [1]. 

In this Letter we propose exact probabilities for a two- 
dimensional Ising system to reach this metastable stripe 
state. Our argument should apply to any curvature- 
driven coarsening process with a non-conserved scalar 
order parameter, such as the time-dependent Ginzburg- 
Landau equation [1, 0, @. For a square system, the 
probability to get trapped in a vertical or a horizontal 
stripe state was previously estimated [HQ to be close to 
i. Here we exploit the relation between coarsening and 
critical percolation to argue that this probability equals 
i — In 11 = 0.3558 ... for a system with free bound- 
ary conditions; the corresponding freezing probability for 
periodic boundary conditions is found to be 0.3390. . . . 

Our approach is based on two key ingredients: 

(i) Soon after the quench, the characteristic domain 
length scale t becomes substantially larger than the lat- 
tice spacing a while remaining much less than the sys- 
tem size: L, a <^ £ <^ L. At this stage, the domain 
mosaic has reached the critical point of continuum per- 
colation (Fig. [ija)), as previously observed in various 
two-dimensional systems 0, 0] ■ 

(ii) During the coarsening regime, £{t) ^ a, the dynam- 
ics becomes deterministic and domain wall evolution is 
driven only by local curvature [H, [1, |6| . Thus the global 
domain topology does not change once the coarsening 
regime is reached Q. 

These features imply that the ultimate fate of the sys- 
tem is pre-determined once the critical percolation state 
is reached. For instance, if a domain exists that crosses 
the system only horizontally or only vertically, a stripe 
state is reached. Figure [1] shows an example of a verti- 
cal spanning domain shortly after the quench that ulti- 



mately coarsens into a vertical stripe; conversely, a do- 
main mosaic that spans in both horizontally and verti- 
cally coarsens into a ground state. 

While individually the observations (i) and (ii) are 
known, their combined use leads to analytical predictions 
for the fate of two-dimensional Ising ferromagnets. The 
connection with the critical state of continuum percola- 
tion is central as it allows us to use exact results about 
crossing probabilities to determine the final state of the 
Ising system. These quantities are defined as the prob- 
abilities for the existence of a spanning cluster with a 
specified topology. We denote these crossing probabili- 
ties for free and periodic boundary conditions as J- and 
V, respectively. For a critical rectangular system with 
aspect ratio r (ratio of height to width) [l^ , the crossing 
probabilities are non-trivial (ie., strictly between and 



1), universal functions of r 111. Il2l|. Beautiful exact ex- 



pressions for various crossing probabilities were originally 
calculated via conformal field theory 1^ 



14 
17 



Il9| 



and 



some of them have been proved in Refs. 

The connection (i) with critical percolation may seem 
surprising since the initial state (with equal numbers of 
up and down spins) is not critical; on the square lat- 
tice, the up spins percolate when their density exceeds 
Pc ~ 0.5927. Hence initially neither phase percolates. 
However the system quickly approaches the critical state 
of continuum percolation in which the crossing probabil- 
ities are non-trivial. To appreciate this property, con- 
sider for concreteness domain coarsening in the square 
geometry with free boundary conditions. Let JF^ be the 
probability that a spin up domain crosses vertically (with 
no constraint on horizontal crossing); J-'y also equals the 
horizontal crossing probability (because the system is 
square), and further equals the crossing of spin down do- 
mains (due to up/down symmetry). A vertical spanning 
spin up domain exists if and only if a horizontal span- 
ning spin down domain does not. (This last statement 
is predicated on the non existence of singular points in 
which more than two domains meet, see Fig. [2l) Hence 
!Fi, = 1 — JF„, so that = ^, which is indeed non-trivial. 

We begin with the analytically simpler case of free 




FIG. 1: Coarsening in the kinetic Ising model on a 1024 x 1024 square lattice with periodic boundary conditions at times (a) 
200, (b) 1000, (c) 5000, and (d) 2.5 x 10* Monte Carlo steps following a quench from T = cx) to 0. Domains are regions of either 
spin up (gray) or spin down (black) . A spanning spin- up domain, that eventually coarsens into a vertical stripe, is highlighted. 





FIG. 2: (a) A state where 4 Ising domains meet at a single 
point evolves into a state (b) where only two domains can 
meet under zero-temperature single-spin flip dynamics. 



boundary conditions. In this setting, every domain mo- 
saic either spans only vertically, only horizontally, both 
horizontally and vertically (dual spanning), or the mo- 
saic does not contain a spanning cluster. Their respective 
probabilities, Tji^, ^hu, ^hv, and J-hv (the dual-spanning 
and non-probabilities are identical by up/down symme- 
try) therefore satisfy the normalization condition 

Tj^,+Tm + 2Thv^l. (1) 

Moreover, the exact form of J-j^^ is known to be [3, [l^ 



^X.M = ||A3F.(l,l,i^,2;A 



(2) 



where pFq{ai, . . . ap[bn, ... 6^; A) is the generalized hyper- 
geometric function [20|, A = A(r) is defined implicitly by 



A = 



with r 



2K{P) 
K{l-k^) ' 



and K(u) is the complete elliptic integral of the first 
kind 0. 

The corresponding horizontal crossing probability fol- 
lows by symmetry 



(3) 



while the crossing probability for dual spanning satisfies 
^hv{r) - i(l - Tj^^ir) - Tm{r)) = TUl/r). (4) 

The basic relation between Ising domains and critical 
percolation implies that, e.g., the crossing probability 
•^hv coincides with the probability for the Ising system 
to freeze into a vertical stripe state as a function of r. 

To test this basic prediction for the freezing probabil- 
ity, we simulate the kinetic Ising model at zero tempera- 
ture using single-spin flip dynamics with the Metropolis 
acceptance criterion — a spin is flipped if its energy de- 
creases or remains the same as a result of the flip. To 
make this simulation more efficient, a list of "active" 
spins — those whose energy will not increase upon be- 
ing flipped — is maintained and constantly updated dur- 
ing the dynamics. In each update step an active spin 
is chosen at random and flipped. One Monte Carlo step 
corresponds to each active spin flipping once, on average. 

We simulate many quenches from T = oo to on lat- 
tices of dimension (256/r) x 256. For each value of r, 
we performed 2 x 10* simulation runs, with each starting 
from a different random initial condition. We define a do- 
main as a connected cluster of nearest-neighbor aligned 
spins. Clusters are identified using a cluster multilabeling 
method 2]J. To determine whether a quenched system 
ultimately freezes into a stripe state, one should, in prin- 
ciple, simulate until the system ceases to evolve. The fi- 
nal stages of the evolution take a disproportionately large 
amount of CPU time, however, and it is advantageous to 
stop the simulation when the domain mosaic first reaches 
its final state topology. For this Ising system, our simula- 
tions indicate that after 200 Monte Carlo steps, domain 
mosaics have reached the topology of the final state with 
a probability that exceeds 0.998. Thus we may identify 
the state of the system at this early time as the predictor 
of the topology in the final state. 

Figure [3] shows our numerical results for the probabil- 
ity for a specified ultimate fate of an Ising system with 
free boundary conditions for a variety of aspect ratios 
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FIG. 3: Probabilities of various domain topologies in the 
kinetic Ising model following a quench with free boundary 
conditions; vertical stripes (+), dual-spanning configurations 
(□), and horizontal stripes (x). Error bars are about 1/3 
the symbol size. The lattice dimensions are (256/r) x 256 for 
various aspect ratios between and 1. Also shown are the cor- 
responding exact percolation crossing probabilities J-j^^, Thv, 
and 2Thv, respectively, from Eqs. ©-(11}. 



FIG. 4; Probabilities of various domain topologies in the ki- 
netic Ising model following a quench with periodic boundary 
conditions; vertical stripes dual-spanning configurations 
(□), and horizontal stripes (x), with error bars about 1/3 
of the symbol size. The lattice dimensions are (256/r) x 256 
for aspect ratios between and 1. Also shown are the corre- 
sponding exact percolation crossing probabilities Pi,o, Vo,!, 
and 2Px, respectively, from Eqs. ([6j and (O. 



r. These agree well with the exact crossing probabiUties 
that follow from Eqs. ((21)^(111) • For the important special 
case of a square geometry, r — I, Eq. ([2]) can be simplified 
to " 



^7r„(l)=^Ul) = J-^ln I = 0.1779. 



(5) 



from which the probability of the Ising system coarsening 
into a stripe state equals 2T-j^^{\) = 0.3558 . . .. An earher 
numerical estimate for the probability of reaching a stripe 
state [ij is consistent with this exact result. 

For periodic boundary conditions, a parallel set of re- 
sults can be constructed to again connect the ultimate 
fate of the Ising system and percolation crossing probabil- 
ities. The nature of the crossing probabihties is substan- 
tially more complex for systems with periodic boundaries 
because spanning clusters can wrap around the torus 
multiple times in the vertical and horizontal directions. 
There are two types of spanning clusters 
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"Winding" clusters are labeled by their vertical and hor- 
izontal winding numbers, (a, 6). For example, winding 
numbers (0, 1) and (1,0) correspond to a vertical and a 
horizontal stripe, respectively. A spanning cluster that 
wraps around the torus once in the vertical direction and 
once in the horizontal direction can have one of two wind- 
ing number pairs, (1,1) or (1,-1), and gives a diagonal 
stripe configuration when the torus is unwound onto the 
square. The other cluster type is the "cross topology" in 
which a spanning cluster is formed by the union of two or 
more spanning clusters with distinct winding numbers. 

Let Vafiii") denote the crossing probability for a span- 
ning cluster with winding numbers (a, b) to exist on a 
rectangle with aspect ratio r and periodic boundary con- 



ditions. This probability is given by [ij]. 



2^(3i+2)a,(3;+2)fc + -Z^(2/+l)a,(2/+l)fc 



where Zm^n is a shorthand for Zm.n{'^',i")', generally 



(6) 
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^,^2(-g-27rr) 



^ — 7rg{m, /r-\-n r) 



and " nig) = q^^^'^Y[k>i(^ ~ Dedekind 77 func- 

tion [20|. Additionally, the co nfig uration with cross 
topology occurs with probability [14 1 



-Pxir) = - 



(7) 



where Z{gJ;r) = fY.m.nez^f^Jn{g]r). By symme- 
try, Vx (f ) also represents the probability that no cluster 
spans the system. 

To compute the crossing probabilities numerically, it 
is preferable to employ the asymptotic expansions [23| 

VoAr) - l-2p^/4 + p-V2pi2_4p53/4^3^i5^^^ 
27 



Vx{r) 



P 



5/4 



53/6 



where p = e^'^/^'', rather than evaluating the special 
functions directly. These expansions provide an excel- 
lent approximation for the entire range of aspect ratio 
< r < 1 [231]. 
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Again, a domain mosaic characterized by winding 
numbers (0,1) [or (1,0)] occurs with probabihty Po,i (or 
Pifi), as given by Eq. ([5]), and coarsens into a vertical (or 
a horizontal) stripe state. Similarly, a mosaic with cross 
topology occurs with probability 2Vx given by Eq. ([7]), 
and coarsens directly into the ground state (either all 
spins pointing up or all pointing down). However, a do- 
main mosaic can also reach the ground state by the in- 
direct route of first forming a diagonal stripe state with 
non-zero winding numbers in both directions. As found 
previously for the specific case of the (1,1) stripe, such 
states are long lived [l|; namely, they reach the ground 
state at a time scale that is much larger than the typical 
coarsening time 0{L^). 

In Fig. [4] we plot the realizations that evolve to a topol- 
ogy with winding number (0, 1) or (1, 0), or to the cross 
topology for a variety of aspect ratios r. These again 
agree well with the exact percolation crossing probabil- 
ities that follow from Eqs. ^ and In the specific 
case of the square system (aspect ratio r — 1), the prob- 
ability of reaching an infinitely along lived stripe state 
is 0.3390 . . .. Because the kinetic Ising model can also 
evolve to diagonal stripe topologies, 7^o,i + ^i,o + 27^ is 
less than 1. 

In conclusion, the probabilities with which Ising ferro- 
magnets freeze into metastable stripe states correspond 
exactly to crossing probabilities in critical percolation. 
Our simulation results for the probabilities to reach a 
specified ultimate fate (Figs.[3]and|4]) agree with theoret- 
ical predictions. Our approach also applies to arbitrarily- 
shaped domains and boundary conditions, and also can 
be used to determine more subtle characteristics, such as 
the distribution of the number of stripes. 

The correspondence with continuum percolation at the 
critical point depends crucially on the initial state con- 
sisting of an equal number of up and down spins. This 
equality always holds if the system is quenched from equi- 
librium at an initial temperature that is above the critical 
temperature Tc- If, however, the system is in a magnetic 
field that is switched to zero simultaneously with quench- 
ing, or if the initial temperature is below Tc, then the 
Ising ferromagnet always freezes into the ground state 
corresponding to the initial majority phase [l]. This be- 
havior immediately follows from the percolation view- 
point. Such a system is above the percolation threshold 
and there always exists a spanning domain of the ma- 
jority phase that percolates in both directions. Conse- 
quently, the system coarsens into the ground state. 

While wc focused on the Ising model with zero temper- 
ature spin-flip dynamics, our theory applies more gener- 
ally to phase ordering kinetics in two-dimensional sys- 
tems with non-conserved scalar order parameter, as well 
as to quenches to sub-critical temperatures. In the latter 



case, the coarsening regime requires that ^ ^ i{t) <C L 
(where the equilibrium correlation length ^ may be ar- 
bitrarily large), and metastable stripe states last a finite 
but large time compared to the coarsening timescale. 
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